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We theoretically discover that unique eigenmodes of skyrmion crystal (SkX) are not only mag- 
netically active to ac magnetic field (H u ) but also electrically active to ac electric field (E u ) in a 
multiferroic chiral- lattice magnet Cu20Se03, which amplifies the dynamical magnetoelectric cou- 
pling between E^ and the spin texture. The resulting intense interference between their electric 
and magnetic activation processes can lead to unprecedentedly large diode effect on the microwave, 
i.e., its absorption by SkX changes up to ~20% when the incident direction is reversed. Our results 
demonstrate that the skyrmion can be a promising building block for microwave devices. 

PACS numbers: 76.50.+g, 75.70.Ak, 75.10.Hk, 75.78.-n 



Skyrmion, a topological vortex-like swirling spin tex- 
ture [1], is now attracting a great deal of interest. 
It was predicted that the skyrmion and its crystal- 
lized form, so-called skyrmion crystal (SkX), are real- 
ized in chiral-lattice magnets without inversion symmetry 
through competition between ferromagnetic (FM) and 
Dzyaloshinskii-Moriya (DM) interactions under magnetic 
field H [H, Quite recently, the SkX phase was in- 
deed observed in metallic B20 alloys such as MnSi 0-Q ? 
Fei-^Co^Si [1, [s[ and FeGe [Toj as well as in an insulat- 
ing magnet Cu20Se03 |lll4l3j by small-angle neutron- 
scattering experiments (SANS) and Lorentz transmission 
electron microscopy (LTEM). 

Since then, several experiments have been performed, 
and have reported intriguing transport properties in 
SkX and electric control of skyrmions with 



spin-polarized current [20|, l2l| in the metallic systems. 



The emergence of spin-driven ferroelectric polarization 
P has been observed in the insulating SkX phase of 
Cu 2 OSe0 3 [H 0, El, and the electric-field control 
of this multiferroic sky rmion texture was experimen- 
tally demonstrated [24[. There the research interest, 
more or less, comes from the possible application to 
next-generation spintronics devices. However a lot of 
researchers are presaging further potentiality in the 
skyrmions. Nevertheless the sorts of experimental works 
are quite limited, i.e., observations by means of LTEM 
or SANS and transport measurements only. 

In this Letter, we theoretically propose a brand new di- 
rection to the research on skyrmions from the viewpoints 
of microwave functionalities and dynamical phenomena 
at GHz frequencies. We discover that collective rota- 
tional and breathing motions of skyrmions in SkX can be 
resonantly activated not only by ac magnetic field (H u ) 
but also by ac electric field (E u ) as unique eigenmodes 
of SkX in multiferroic chiral-lattice magnets. These reso- 
nances amplify the dynamical coupling of underlying spin 
texture with E u and , and the resulting intense in- 



terference between the electric and magnetic activation 
processes can lead to unprecedentedly large directional 
dichroism of electromagnetic (EM) wave in Cu2 0Se03, 
i.e., its absorption by SkX changes up to ~20% depend- 
ing on the sign of its incident direction. This effect is 
enhanced especially at eigenfrequencies of the aforemen- 
tioned skyrmion resonances (~GHz), and can work as 
an efficient microwave diode. Currently, most of the 
microwave- device functions are achieved using designed 
combinations of waveguides, circuits and elements made 
of ferrites with ferrimagnetic order [25J . Our finding pro- 
vides a guideline for designing new microwave devices 
such as magnetically tunable isolator. Our work will be 
a trigger for a broad-based quest for novel functions of 
skyrmions and related spin textures. 

The magnetic structure of Cu20Se03 is composed 
of tetrahedra of four Cu 2+ (S= 1/2) ions as shown in 
Fig. [T]fa). Recent powder neutron diffraction [26] and 
NMR [27[ experiments suggested that three-up and one- 
down (uuud) type collinear spin arrangement is realized 
on each tetrahedron below T c ~58 K. We regard this four- 
spin assembly as a magnetic unit, and treat it as a clas- 
sical vector spin rrti whose norm m is unity. We employ 
a classical Heisenberg model on a cubic lattice |28|-[30| to 
describe the magnetism in thin specimen of Cu2 0Se03, 
which contains the FM-exchange interaction, the DM 
interaction [3l| and the Zeeman coupling to the external 
H normal to the plane. The Hamiltonian is given by, 



Ho = -j ^2 m i • m j - p 

<i,j> i,7 



m , x rn 



(1) 



where g=2, and 7 runs over x and y. We set the ratio 
.D/J=0.09, for which the periodicity in the SkX phase 
becomes ^99 sites. Since the distance between adjacent 
Cu-ion tetrahedra is ~5 A, this periodicity corresponds 
to the skyrmion diameter of ~50 nm in agreement with 
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FIG. 1: (color), (a) Spin structure of Cu2 0Se03, composed 
of tetrahedra of four Cu 2+ -ions (£=1/2) with three-up and 
one-down spins, (b) Phase diagram of the spin model (pQ). 
Here HL, SkX, and FM denote helical, skyrmion-crystal, and 
ferromagnetic phases, respectively, (c) Spin structure in the 
SkX phase, which possesses a six-fold rotation axis, 6, and six 
two- fold rotation axes followed by time reversal, 2 / . Arrows 
represent the in-plane spin components, (d) Under if || [110], 
the system becomes polar along [001] and the emergence of 
ferroelectric polarization P||[001] is allowed, (e) Symmetry 
axes in the Cu2 0Se03 crystal, which belongs to the P2i3 
space group: three- fold rotation axes, 3, along (111) and two- 
fold screw axes, 2i, along (100). (f) Real-space configura- 
tions of local electric polarizations pi in the skyrmion under 
ff||[110]. 



the LTEM observation [ll|. All the spin textures con- 



sidered here are slowly varying, which can be described 
by a continuum spin model. It justifies our treatment 
based on a lattice spin model after dividing the space 
into square meshes and the coarse graining of magneti- 
zations. 

We first analyze the model (pQ) using the replica- 
exchange Monte-Carlo technique, and obtain phase di- 
agrams at low temperature (T) shown in Fig. [Tfb). 
The SkX phase emerges in the range 1.875 x 10 -3 < 
\gfJLB^oH z /J\ < 6.3 x 10 -3 , sandwiched by the helical 
and FM phases in agreement with the experiment for 



thin-plate samples [llj. The skyrmions are crystallized 



into a triangular lattice and the magnetic moments mi 
direct antiparallel (parallel) to if at the center (periph- 
ery) of each skyrmion as shown in Fig. [U^c). 

For the SkX state formed under if ||[110] as shown in 
Fig.QJd), the emergence of P||[001] perpendicular to the 
net magnetization M(\\H) is expected from symmetry 
consideration. As shown in Fig. Die), the crystal struc- 



ture of Cu20Se03, which belongs to a non-polar space 
group, P2i3, possesses three-fold rotation axes, 3, along 
(111), and 2i-screw axes along (100). The spin texture 
in the SkX phase is also non-polar with a six-fold rota- 
tion axis, 6, along if, and two- fold rotation axes followed 
by time reversal, 2', normal to if as shown in Fig. [U^c). 
When the SkX sets in under if ||[110] on the Cu 2 OSe0 3 
crystal, most of the symmetries should be broken, and 
only the 2^-axis ( 1 1 [001] ) normal to if survives as shown 
in Fig. QJd). Consequently, the system becomes polar 
along [001]. Indeed the emergence of P||[001] under 
if || [110] was experimentally observed [23]. Here we de- 
fine the Cartesian coordinates, a? || [110], 2/||P|| [001], and 
z\\M || [110] shown in Fig. DJd), for convenience of the 
following formulations. 

The net magnetization M and the ferroelectric polar- 
ization P are given by sums of the local contributions 
as M=j^^^ 1 m i and P=^y Ym=i Ph respectively, 
where the index i runs over the Cu-ion tetrahedra with 
uuud spin pair, N is the number of the tetrahedra, and 
V(=1.76x 10 -28 m 3 ) is the volume per tetrahedron. Be- 
cause of the cubic symmetry, the local polarization pi 
from the ith tetrahedron is given using the spin compo- 
nents mj fl , m^, and m^ c in the P2i3 setting as, 



Pi 



(Pia,Pib,Pic) = A (m ib m ic , m ic m ia , m ia m ib ) . (2) 



We can easily evaluate the local contributions pi and rrti 
from each tetrahedron in the ferrimagnetic phase where 
all the tetrahedra give uniform contributions. Then the 
coupling constant A is evaluated as A=5.64 x 10 -27 /iCm 
from the experimentally measured -P[ooi]— 16 /iC/m 2 in 
the ferrimagnetic phase under if 1 1 [111] at 5 K [ll| . 

Because of this strong coupling between magnetism 
and electricity, collective oscillations of this SkX can be 
activated not only magnetically by ac magnetic field if 00 
but also electrically by ac electric field E u . As demon- 
strated below, with the special configuration of P±M , 
both the if u and E u components of EM wave propagat- 
ing along PxM can activate common oscillation modes. 
To see this, we calculate dynamical magnetic and dielec- 
tric susceptibilities, 



lv± a 

PqH u q 



(3) 



by numerically solving the Landau-Lifshitz- Gilbert equa- 
tion using the fourth-order Runge-Kutta method. The 
equation is given by 



drrti 
dt 



x iff* + — m, x — — - 

rn dt 



(4) 



where aQ(=0.04) is the Gilbert-damping coefficient. The 
effective field iff ff is calculated from the Hamiltonian 
H=Uo+H'(t) as Hf = -dH/ drrti. Here the first term 
Ho is the model Hamiltonian (pQ), while the perturba- 
tion term H'(t) represents a short rectangular pulse of 
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FIG. 2: (color online), (a) Imaginary parts of the dynamical 
magnetic susceptibilities, Imx™ and Imx™ m , as functions 
of frequency uo at g/^B/^oHz/ J=3-75x 10 -3 . In Ref. [33], the 
strong (weak) resonance in Imx™ was ascribed to the coun- 
terclockwise (clockwise) rotation mode, while the resonance 
in Imx™ was to the breathing mode, (b) Imaginary parts 
of the dynamical dielectric susceptibilities, ImxTz and Imx^,, 
as functions of uj. Above mentioned three magnetically active 
resonances can be seen in the spectra of dielectric suscepti- 
bilities at the same frequencies, indicating their simultaneous 
electric activities. 



magnetic field or electric field. After applying the pulse 
at £=0, we calculate M(t) and P(t), and obtain their 
Fourier transforms M£ and P£. The calculations are 
performed using a system of A/"=288x288 sites with the 
periodic boundary condition. 

In Fig [2ja), we display imaginary parts of the cal- 
culated dynamical magnetic susceptibilities, Imx™(o;) 
and Imx£TM, for H u \\y and H u \\z. We also plot 
imaginary parts of the calculated dielectric susceptibil- 
ities, lmx e z e z (oo) and lm.Xy e y (co>), for E u \\z and E u \\y in 
Fig[2jb). In Imx™, we find a strong resonance active to 
H u \\y at o; R /J=6.12xlO- 3 , which was ascribed to the 
counterclockwise rotation mode where all the skyrmion 
cores in the SkX uniformly rotate in the counterclock- 
wise fashion [32|, [33[. This rotation mode can be seen 
also in the spectrum of lvax e zz as a peak at the same 
frequency indicating its simultaneous electric activity to 
E^\\z. The spectrum of Imx™ has one more resonance 
at higher frequency cjr=1.135 x 10 _2 J, which was as- 
cribed to another rotation mode with opposite rotational 
sense, i.e., the clockwise rotation mode [32|, [33j. We can 
see a very tiny peak in Imx e z e z at the corresponding fre- 
quency, indicating its weak electric activity. On the other 
hand, the spectrum of Imx™ has a single resonance ac- 
tive to H u \\z at cjr/J=7.76x10 -3 , which was ascribed 
to the breathing mode where areas of all the skyrmions 
in the SkX oscillatory expand and shrink in a uniform 
way [32|, [33| . Again this mode is simultaneously active 
to E^Wy, and the corresponding peak can be seen in 
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FIG. 3: (color online). Configurations of the microwave 
and E u components, for which the nonreciprocal directional 
dichroism is expected when P\\y and M\\z with PA_M, and 
K w ||PxM(||aj); (a) E u \\z and H u \\y, and (b) E u \\y and 
H u \\z. Cartesian coordinates aj||[110], y||[001], and z\\[110] 
are defined as shown in the figures where the z axis is parallel 
to H. 



Imx^ at the same frequency. Recent microwave exper- 
iment found clear absorptions at these spin-wave reso- 
nances, while absorptions at off-resonant frequencies turn 
out negligibly small [34[. 

The presence of collective modes active to both E u 
and H u is nothing but the source of interesting mi- 
crowave activity. From Maxwell's equations, we can de- 
rive the relation H^WK^xE" for the EM wave. This 
relation indicates that relative directions of H u and E" 
are determined by the propagation vector K u , and their 
relationship should be reversed upon the sign reversal 
of K u . When the lineally polarized EM wave with 
E^\\z and H^Wy propagates parallel (antiparallel) to 
(P x M)||as as shown in Fig.[3ja), where sgn[Rei ; C u; ]=+l 
(sgn[ReK u; ]=-l) with K UJ =K UJ x, the oscillation of P 
induced by E u and that of M by H u contributes in a 
subtractive (an additive) way to the collective oscillation, 
which results in weaker (stronger) absorption of the EM 
wave. Such a nonreciprocal absorption of EM- wave is ex- 
pected also for E u \\y and H u \\z as shown in Fig. E9(b). 

To study the microwave absorption and NDD quan- 
titatively, we start with the following Fourier-formed 
Maxwell's equations for materials with M and P [35[, 



uB" = K u x E", -cuD" = K^ x H u 



(5) 
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FIG. 4: (color online). Calculated absorption coefficients, 
a+(cj) and ol—{uS) for microwaves with sgn(ReX u; )=H-l and 
sgn(ReK UJ ) = — 1, respectively, at several values of H z in the 
cases of (a) H u \\y and E u \\z and (b) H u \\z and E u \\y. 



where 

B" = f jLo(fi, 00 H u + M u '), D u = eoe^E" + P w ,(6) 

r eo~ 



V Mo 

P" = e r e ^)^+X em (^)v / ^o"^ a; . 



(7) 

(8) 



As discussed above, the NDD of microwaves with 
K UJ =K UJ x is expected for the following configurations 
of H u and E u \ (i) H u \\y and and (ii) H^Wz 

and E/^Hy. We introduce the complex refractive in- 
dex N (uo)=n(uo) + m(cj) which is related with K u as 
K u =^N((j). We solve Eqs. (j5j) and obtain, 



7V(cj) 



« 00 + xSM]k 00 + C( w )l 

^n(Re^)[x^ e H+x e ^H]/2, 



(9) 



and 



~ ^/[e»» 0O + x^H]lM« 0O + xSrH] 

+ Sff n(Re^)[x^(c) + x e y TH]/2, (10) 

for Cases (i) and (ii), respectively. These expressions con- 
tain the sign of ReK" indicating direction dependence of 



the microwave absorption because the absorption coeffi- 
cient a{uj) is related to N(uo) as, 

afcj) = ex cjImVV(cj), 

c 

and the absorption intensity is given by 
I(u)=Iq exp[— a{uj)l] where I is the sample thickness. 
The nonreciprocal absorption Aa(uj) = &+(uj) — a_(cj) 
represents the magnitude of NDD where a + and a_ are 
absorption coefficients for microwaves propagating in 
the positive and negative directions, respectively. 

In order to evaluate N(lj) and a±(u) quantitatively, 
we need to calculate not only the dielectric and magnetic 
susceptibilities but also the following dynamical ME sus- 
ceptibilities, 



x e Q » 




(11) 



For the values of and in Eqs. (|9|) and ([T0]h we 
assume an isotropic dielectric tensor, i.e, e^=e^=e°° 
for simplicity, and set e°°=8 according to the dielectric- 
measurement data 22|, 36[. In turn, we take /i^=/i^=l 
for permeability. The value of J is set to be J=l meV 
so as to reproduce the experimental T c for the SkX- 
paramagnetic transition. 

In Fig. Ufa) [Fig. Efb)], we display calculated lj 
dependence of the absorption coefficients, &+(uj) and 
a -(w), at several values of H z for E u \\z and H u \\y 
[E u \\y and H"\\z\. The out-of-plane E u \\z and in- 
plane H^Wv activate two rotation modes with oppo- 
site senses [32|, [33[, i.e., lower- lying counterclockwise and 
higher- lying clockwise modes, which give two spectral 
peaks in Fig. 21(a). We find that for the lower- lying reso- 
nance, the Aa(uj) increases as H z increases, and reaches 
more than 0.25 cm -1 , which corresponds to a relative 
change Aa/a ave =2(a + — a_)/(a + -f a_)~20 % at max- 
imum. On the other hand, both the in-plane E^Wy and 
out-of-plane H u \\z activate a breathing mode [32I I33I]. 
which give a single spetral peak in Fig. [U(b). Again the 
Aa(uj) increases with increasing H z , and reaches approx- 
imately 0.14 cm -1 , corresponding to Aa/a ave of 10 %. 
These values do not depend on the value of ac- Such 
a huge directional dichroism is quite rare for any fre- 
quency range [HI, [3tM3Q| ] . and has never been realized 
at GHz frequencies. This is because most of the multi- 
ferroics with simple magnetic orders have resonant fre- 
quencies much larger than microwave frequencies due to 
the large spin gaps. In turn, the long-period skyrmion 
textures have small spin gaps and their nontrivial col- 
lective modes with GHz frequencies enable us to achieve 
interesting microwave functions. 

In summary, we have theoretically predicted that 
the SkX phase of chiral-lattice insulator Cu2 0Se03 
shows the enhanced diode effect on linearly polarized 
microwaves as a consequence of interference between 
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magnetic and electric responses from the multiferroic 
skyrmion texture. Our prediction demonstrates that the 
skyrmion and the SkX host interesting dynamical phe- 
nomena in the microwave frequency regimes in addition 
to the peculiar transport and spintronics phenomena. 
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